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Abstract We investigate the cold nuclear matter(CNM) ef- 
fects on dijet productions in high-energy nuclear collisions 
at LHC with the next-to-leading order perturbative QCD. 
The nuclear modifications for dijet angular distributions, di- 
jet invariant mass spectra, dijet transverse momentum spec- 
tra and dijet momentum imbalance due to CNM effects are 
calculated by incorporating EPS, EKS, HKN and DS param- 
etrization sets of parton distributions in nucleus . It is found 
that dijet angular distributions and dijet momentum imbal- 
ance are insensitive to the initial-state CNM effects and thus 
provide optimal tools to study the final-state hot QGP ef- 
fects such as jet quenching. On the other hand, the invari- 
ant mass spectra and the transverse momentum spectra of 
dijet are generally enhanced in a wide region of the invari- 
ant mass or transverse momentum due to CNM effects with 
a feature opposite to the expected suppression because of 
the final-state paiton energy loss effect in the QGP. The dif- 
ference of EPS, EKS, HKN and DS parametrization sets of 
nuclear parton distribution functions is appreciable for dijet 
invariant mass spectra and transverse momentum spectra at 
pH-Pb collisions, and becomes more pronounced for those at 
PbH-Pb reactions. 

1 Introduction 

In relativistic heavy-ion collisions with a large amount of 
energies deposited in the collision center quarks and glu- 
ons confined in nucleons should be liberated to form a new 
kind of matter - quark gluon plasma (QGP). Hard processes 
with large momentum transfer have long been regarded as 
good probes of the QGP created in high-energy nuclear re- 
actions [1] and an important probe among them is the jet 
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quenching, or the energy loss of an energetic parton in QCD 
matter [2-5]. With the upgrading facilities at RHIC and es- 
pecially the unprecedented colliding energies at LHC, mea- 
surements of leading particle production at large transverse 
momentum have been extended to reconstructed jet produc- 
tion, and several theoretical studies of jet production in high- 
energy nuclear collisions have been made by investigating 
hadron multiplicity inside a jet [6], jet shape [7, 8], inclu- 
sive jet cross section at the next-to-leading order (NLO) [9], 
and Z°/f -tagged jets [10]. Recently ATLAS and CMS have 
made the first successful measurement of jet production in 
Ah-A reactions at LHC by observing dijet asymmetry in Lead- 
Lead collisions [11, 12], which has aroused intense theo- 
retical investigations [13-16]. It has been demonstrated that 
jet quenching, or final-state parton energy loss in the QGP 
should be responsible for the large dijet momentum imbal- 
ance in nuclear reactions at LHC [17] with some caution of 
the background fluctuation effect in heavy-ion collisions [18]. 

In heavy-ion collisions because the initial-state colliding 
objects are nuclei instead of nucleons, cold nuclear matter 
(CNM) effects such as shadowing effects, anti-shadowing 
effects, EMC effects etc. will manifest themselves in hard 
processes, and it has been shown that CNM effects will play 
an important role in productions of leading hadrons [19], di- 
rect photons [20-22], Drell-Yan dilepton [19, 23, 24], heavy 
flavor mesons [25-27] as well as inclusive jets [9, 28] and 
other hard processes in nuclear collisions. It will be of great 
interest to see how large CNM effects will be in dijet produc- 
tion of high-energy nucleus collisions, which will then help 
us understand the final-state jet quenching effect more pre- 
cisely. It has been known that for hard processes in nuclear 
collisions with the assumption of the leading twist factor- 
ization of perturbative QCD (pQCD), the main CNM effects 
can be phenomenologically parameterized the difference be- 
tween parton distribution functions (PDFs) in nucleon and 
nuclear parton distribution functions (nPDFs) [29]. Other 
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theoretical models have also been developed to pinpoint the 
underlining mechanisms of different CNM effects with many- 
body QCD and related phenomenologies [30, 31]. In this 
paper we will investigate CNM effects on dijet production 
in p+Pb and Pb+Pb collisions at LHC at the NLO perturba- 
tive QCD by incorporating 4 parametrization sets of nPDFs 
- EPS [29], EKS [32], HKN [33] and DS [34]. Four di- 
jet observables — dijet angular distributions, dijet invariant 
mass spectra, dijet transverse momentum spectra and di- 
jet momentum imbalance with CNM effects in pH-Pb and 
PbH-Pb collisions at LHC will be calculated as well as the 
nuclear modification factors for these dijet observables due 
to CNM effects. 

Our work is organized as follows: in Section 2, we give 
the formalism of dijet productions in nucleon-nucleon col- 
lisions with perturbative QCD. Dijet observables for pH-Pb 
and PbH-Pb colhsions with CNM effects at LHC by utiUz- 
ing 4 different sets of nPDFs and their physics implications 
are studied in Section 3. At last we give a brief summary in 
Section 4. 



2 Dijet Production in nucleon-nucleon Collisions 

Dijet productions in elementary nucleon-nucleon collisions 
provide the base line for dijet productions in relativistic heavy- 
ion collisions. With the perturbative QCD factorization ap- 
proach [35, 36], we can obtain leading-order (LO) invariant 
transverse momentum spectrum for dijet in elementary col- 
lisions 
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where fa,b{xa.b) is the parton distribution function in nu- 
cleon, and Xa, are the momentum fractions carried by 
participating partons from nucleons. They have relationships 
with transverse energy Et and rapidity yi, y2 of final parti- 
cles as 



(2) 



da /dt is the elementary scattering cross section for partons 
at tree level. And t is the Mandelstam variable. 

At next-to-leading order (NLO), jets can be defined by 
jet-finding algorithms with a radius parameter in rapidity y 
and azimuthal angle plane, R = \J Ay'^ +40^ [37] xhe 
cross section of dijet production at NLO can be expressed 
as [38], 
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Fig. 1 The angular distributions of dijet are evaluated at NLO and 
compared to experimental results fonn DO collaboration in p+/5 colli- 
sions at ,/s = 1 .96 TeV [50] and CMS collaboration in p+p collisions 
atv^ = 7 TeV [51]. 



Here, Vj represents the physical quantity of final state, and 
, Eji, yi, 0, are the four-momentum, transverse energy, 
rapidity, azimuthal angle of the /— th(/ — 1, 2, 3) particle, 
respectively. The first term on the right-hand denotes the 
contribution from 2 — >^ 2 processes including NLO virtual 
corrections. The second term is the contribution from 2^3 
processes. And the functions S contain the jet finding algo- 
rithm. In this paper we will estimate Vj as angular, invariant 
mass and final transverse momentum of dijet by utilizing the 
EKS framework of a NLO calculation of jet productions in 
hadron-hadron collisions [37-39]. 

The rapidity of a dijet system can be defined as 



yi -\-y2 * bi - 



yi] 



(4) 



and the angular of system is related to rapidity as j = e 
For massless partons, x can be written as 



2y* 



x = 



COS 9* 
cos 6* 



(5) 



where 9* is the polar scattering angle of the outgoing jets in 
the dijet center-of-mass frame [40]. 
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Fig. 2 Compaiison between dijet invariant mass spectrum at the NLO 
with the jet size R = 0.4 and ATLAS data in p+p collisions at y'i = 
7 TeV [53]. 



Fig. 3 Comparison between dijet asymmetry distribution at the NLO 
with the jet size R = 0.5 and the CMS measurement in p+p collisions 
at v^ = 7TeV [12] . 



We can also define the dijet invariant mass Mjj as the 
invariant mass [(L/'n )^]'''^ of all particles in the two jets. 
At LO, it has a form 



Mjj = 2E^[l+cosh{yi^y2)]. 



(6) 



The invariant mass cross section of dijet as well as dijet an- 
gular distribution has provided a good tool to test predictions 
of perturbative QCD and predictions beyond the Standard 
Model(SM) such as quark compositeness, extra spatial di- 
mensions [41^8]. 

Our results for angular distributions at NLO shown in 
Figure 1 are compared to data from DO collaboration and 
CMS collaboration with the utilizing of the CTEQ6.1M par- 
ton distribution functions (PDFs) [49]. DO detector in Run 
II of the Fermilab Tevatron Collider measured the angular 
distribution with the jet size R = 0.7, \yjj\ < 1 at y/s — 
1.96 TeV [50]. And the CMS collaboration measured that 
with the jet size R = 0.5, \yjj\ < 1.11 at ^ = 7 TeV [51]. 
We confront the theoretical results in p+p collisions with DO 
measurements in the invariant mass intervals of 250 GeV< 
Mjj< 300 GeV , 600 GeV< Mjj < 700 GeV, 900 GeV< 
Mjj < 1000 GeV, and the numerical simulations in pH-p reac- 
tions with CMS in the invariant mass intervals of 500 GeV < 
Myy < 650 GeV, 1100 GeV<Myy < 1400 GeV, 1800 GeV< 
Mjj < 2200 GeV. We observe a good agreement between 
theory and experimental data in different intervals of in- 
variant mass. The most symmetrical production at yi — y2, 
namely = 1 is the largest. When the scattering angle 0* 
is very small, the angular distribution is proportional to the 
Rutherford cross section [52]. Hereby, the dijet angular dis- 
tribution approaches a constant at large X- 



The ATLAS collaboration at the LHC has published re- 
cent results about the invariant mass spectrum of dijet in 
the bin of maximum rapidity of the two leading jets, \y\miix 
=max(|3'i|, 1^2 1) in pH-p collisions at y/s = 7 TeV [53]. We 
test the dijet invariant mass cross section in the rapidity bin 
of 0.3 < \y\max < 0.8 with the jet size R = 0.4 in Figure 2. 
The NLO pQCD theory can describe the differential yield 
well in the mass range from 200 GeV to 650 GeV where the 
cross section falls by more than two orders of magnitude . 

At the NLO the exact transverse momentum balance be- 
tween Eti and £7-2 should be broken because of additional 
gluon radiation in 2 3 processes. To quantify the trans- 
verse momentum imbalance between the leading jet and the 
subleading jet one defines the dijet transverse momentum 
asymmetry A j as. 



£7-1 —Et2 
Ej 1 +Et2 



(7) 



The A J distribution can be expressed from the dijet trans- 
verse energy spectrum as follows [16]: 

da /-A lniiix «'2max rEr I lEj 2 



da[ET\{Aj,ET2)] 
dy\dy2dEj idEj2 



(8) 



where /sri =£72(1 —-^j) is the transverse energy 

for leading jet. At the LO we always have Aj ~Q due to the 
equal transverse energies of two jets. 

The CMS collaboration has carried out the asymmetry 
distribution of dijet in pH-p collisions at — 7 TeV [12], 
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requiring the leading jet with the transverse energy of £7-1 > 
120 GeV, and the subleading jet with £■7-2 > 50 GeV in op- 
posite Acj) — 101 — 02| > 27i:/3. The numerical result from 
Eq. (8) is given in Figure 3 and confronted with the exper- 
imental data from CMS. The result from the NLO pQCD 
describes the CMS dijet asymmetry data very well. 



3 Dijet production in heavy-ion collisions with CNM 
effects 

Extrapolating the cross sections in pH-p collisions to p+A and 
A+A, one can replace /^/^ with different sets of nPDFs ( 
EPS09 [29], EKS98 [32], HKN [33], DS [34] ). They mod- 
ify the parton distribution in nucleon with shadowing (anti- 
shadowing) factors which reflect the modification for prob- 
ability of initial participating partons found in the nucleus. 
At the same time, isospin symmetry for bound protons and 
neutrons, EMC effect and Fermi motion in nucleus are also 
included. For a nucleus A with Z protons, the factorization 
form is 

fa/Aix,Q^) = Ra/A{x,Q^)[jfa/p{x,Q^) 



+ (1-t)/«/«(^,G')] 



(9) 



To make the NLO computations of dijet in relativistic heavy- 
ion collisions we combine the NLO calculation of dijets 
in hadron-hadron collisions in Eq. (3) with the nPDFs of 
EPS09 NLO, HKN NLO and DS NLO sets, respectively. 
Because EKS98 does not have NLO nPDFs the LO parametriza- 
tion of EKS98 nPDFs will be used when calculations with 
EKS are carried out. 

We also need insert thickness function Ia (b) in the whole 
reaction plane, which indicates the number of participating 
nucleons for unit area in the nucleus . With the Glauber 
model the spectrum for the dijet production in p+A colli- 
sions is given by 
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and for A + A colhsions, the spectrum is 
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In order to investigate the CNM effects on angular dis- 
tributions in p+A and A+A collisions intuitively, we define 
the nuclear modification factors of angular distributions. 
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Fig. 4 Dijet angular distributions at NLO are extended to p+Pb colli- 
sions at ,/s = 8.8 TeV (top) and Pb+Pb collisions at y'i = 2.76 TeV 
(bottom) with the cone size R = 0.4. Different sets of nPDFs given by 
EPS, EKS, HKN, DS are considered. 



taking the angular distribution in pH-p collisions as a base- 
Une. 

Similarly, the nuclear modification factors for invariant 
mass spectrum d^ GpA^AA) / dMjjd\y\,„a_-(, final transverse mo- 
mentum spectrum d^GpAi^AA) /dETidEri, and momentum im- 
balance dOpAi^AA) /'^pA{AA)dA J can be written respectively as 
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Fig. 5 Nuclear modification factors of dijet angular distributions for 
p+Pb collisions at y(s = 8.8 TeV (top) and Pb+Pb collisions at </s = 
2.76 TeV (bottom) with different nPDFs. 



Fig. 6 Dijet invariant mass spectra at NLO for p+Pb collisions at^/s = 
8.8 TeV (top) and Pb+Pb collisions at y/s = 2.76 TeV(bottom) in the 
rapidity bin 0.3< \y\„u,x <0.8 with different nPDFs. 



pAy _ dapA(AA)/<^pA{AA)dAj 

%AiAA)- da„,/ap,dAj ■ ^'"^ 

The parameter {Nbi„ary} in the denominator represents pairs 
of participating collisions. In p+A collisions it is decided by 

{Nbincry) = jd^htAih). In A+A collisions it is {Nhi„ary) = 

Now we predict dijet productions in minimum bias p+Pb 
and Pb+Pb collisions at LHC due to the CNM effects. Please 
notice that for dijet productions in Pb+Pb collisions the final- 
state hot QGP effects such as jet quenching in the hot QCD 
medium are not included in the investigations presented in 



this paper, and in the following we focus only on the CNM 
effects on dijet productions in heavy-ion reactions. 

Firstly we study the angular distributions in the mass 
interval of 600 GeV to 700 GeV and < 1 for p+Pb 
collisions at ^/s — 8.8 TeV, for Pb+Pb collisions at = 
2.76 TeV. The jet radius is set to 0.4. Both angular distri- 
butions in p+Pb and Pb+Pb shown in Figure 4 have iden- 
tical trends with that in p+p collisions. Results with differ- 
ent nPDF parameterizations EPS, EKS, HKN and DS, show 
very small deviations. 

Figure 5 shows nuclear modifications of dijet angular 
distributions with different sets of nPDFs in p+Pb and Pb+Pb 



6 



1.4 
1.2 
1 

^<0.8 
r 

0.6 
0.4 
0.2 




\fs=8.8TeV p+Pb 
0.3<|y| <0.8 R=0.4 

max 



« EPS 
• EKS 

' HKN 
■ DS 

tSSi EPS09 error band 



200 300 400 500 600 700 800 
Mjj(GeV) 



1.4 
1.2 
l" 

^=^<0.8 

DC 

0.6 
0.4 
0.2 




- 

- 




_ Njs=8.8TeV p+Pb 




0.3<|y| <0.8 




_ '■' 'max ^ 


R=0.2 


• 


R=0.4 




R=0.8 




< 1 . . < < 1 < < < 



200 300 400 500 600 700 800 
Mjj{GeV) 




s=^S0.8t 

DC 

0.6 

0.4; 
0.2; 

0^ 



^=2.76 TeV Pb+Pb 
0.3<|y| <0.8 R=0.4 

' max 



X EPS 
« EKS 
' HKN 
DS 

tJJJS EPS09 error band 



200 ' ' 300 ' ' 400 ' ' 500 ' ' 600 ' ' 700 ' ' 800 



DC 



1.4 
1.2, 
1 

0.8 
0.6 
0.4 
0.2 








' • • 1 1 


* * » i I 1 . . 




, 


- ^=2.76 TeV 


Pb+Pb 


0.3<|y| <0.8 

max 


" R=0.2 




R=0.4 




' R=0.8 


~, < 1 , , < < 1 < < < , 1 < 





MJGeV) 



200 300 400 500 600 700 800 
Mjj(GeV) 



Fig. 7 Nuclear modification factors of dijet invariant mass spectra for 
p+Pb collisions (top) and Pb+Pb collisions (bottom) with different 
nPDFs. 



Fig. 8 Nuclear modification factors of dijet invariant mass spectra for 
different jet cone sizes in a fixed rapidity bin 0.3< \y\max <0.8. Top 
plane is the modification for p+Pb collisions and the bottom is for 
Pb+Pb collisions. 



collisions at LHC. It reveals that the nuclear modifications 
are independent of the scattering angular x and invisible 
over the whole range for both p+Pb and Pb+Pb collisions. 
This interesting feature results from the constrained kine- 
matics regions probed by dijet angular distributions for fixed 
Mjj. At LO from Eq. (2), Eq. (4) and Eq. (6), we can obtain 
Xa*xt, = Mjj/s, which imposes a constraint for the momen- 
tum fractions and xi,. In particular when yjj ^ 0, at LO 
we find Xa ^ X[, ^ Mjjj ^s, which means that the momen- 
tum factions in nPDF show very weak dependence on the 
dijet angular %. We note that in A+A reactions the insensi- 
bility of dijet angular distribution to the initial-state CNM 



effects implies the dijet angular distribution is a good tool 
to probe the final-state hot QCD matter effects such as jet 
quenching in the QGP. 

The invariant mass spectra in a rapidity bin 0.3 < \y\mux < 
0.8, with jet size R = 0.4 for p+Pb colHsions at ^ = 8.8 TeV, 
and Pb+Pb collisions at y/s = 2.76 TeV are plotted in Fig- 
ure 6. The cross sections with four sets of nPDFs fall con- 
sistently with the dijet mass. 

Shown in Figure 7 are nuclear modification factors of 
dijet invariant mass spectra with 4 different sets of nPDFs. 
The overall trend is that the nuclear modification factors 
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Fig. 9 Nuclear modification factors of dijet invariant mass spectra for 
different rapidity bins with a fixed cone size R=OA. Top plane is the 
modification for p+Pb collisions and the bottom is for Pb+Pb colli- 
sions. 



are greater than 1 due to the anti-shadowing effect in the 
kinematical region of Mjj, except that results with EPS are 
slightly suppressed at very large Mjj in Pb+Pb reactions. 
The top plane is the modification factor in p+Pb collisions at 
y/s = 8.8 TeV, and the CNM effects are not very visible. The 
bottom plane is the modification factor in Pb+Pb collisions 
at y(s = 2.76 TeV. In the mid-rapidity region, we can roughly 
estimate atLOM/y ^ lEj viaEq. (6), andEj from 
Eq. (2). Taking HKN set as an example, when Mjj lying 
above - 150 GeV at = 8.8 TeV, and above - 50 GeV at 
y/s — 2.76 TeV the momentum fractions of the initial-state 



partons are in the anti-shadowing region. This leads the en- 
hancement in the range shown in Figure 7. There are devi- 
ations between different sets of nPDFs for both p+Pb and 
Pb+Pb collisions. Therefore the invariant mass spectrum of 
dijet may provide a convenient physical quantity to distin- 
guish different parametrizations of nPDFs. Note that EPS09 
provides the tool to study the theoretical uncertainty of the 
nPDFs parametrization [29] and we calculate the resulting 
error band of R^g on dijet invariant mass distribution as il- 
lustrated in Figure 7 by utilizing 30 different parametrization 
sets of NLO nPDFs from EPS09. One can observe that with 
the error band the deviations between results of EKS, EPS 
and DS are rather small whereas that of HKN gives a con- 
siderably larger enhancement and shows a visible distinction 
with those of EKS, EPS and DS. 

We further investigate the dependence of nuclear mod- 
ification for invariant mass spectrum on the jet radius and 
rapidity bin with EPS nPDF 

Figure 8 demonstrates the sensitivity of the modification 
factor for dijet invariant mass spectrum to the jet size. We 
check three jet cone sizes R — 0.2, 0.4, 0.8 in a fixed rapidity 
bin 0.3 < \y\„uix < 0.8. The modification factor in p+Pb col- 
lisions is almost independent on the jet cone radius R, while 
the modification factor in Pb+Pb collisions demonstrates a 
slight dependence on the jet radius. This phenomenon be- 
cause of the CNM effects is opposite to the nuclear modi- 
fications for single jets [8] and Z"//* -tagged jets [10] with 
final-state hot-dense medium effect in Pb+Pb collisions, where 
one observes strong dependence of R'^ on the jet radius 
R and jet quenching effect decreases with larger jet cone 
size because more energy carried by radiated gluon will fall 
within the jet area when the jet radius R becomes larger. 

Figure 9 illustrates the sensitivity of the modification 
factor for dijet invariant mass spectrum to the rapidity bin. 
We vary rapidity in four bins of \y\„iax < 0.3, 0.8 < \y\max < 
1.2, 1.2 < \y\mux < 2.1, 2.1 < \y\max < 2.8 with the jet cone 
size R — 0.4. The sensitivity of nuclear modification in p+Pb 
collisions to the different rapidity bins is rather modest, and 
we find that the modifications in central rapidity bins are 
stronger than that in 2.1 < \y\„,ax <2.8. This kind of depen- 
dence on rapidity region in Pb+Pb collisions becomes ap- 
parent. The closer the rapidity bin is to central region, the 
stronger the nuclear modification is, and for the largest ra- 
pidity region 2.1 < \y\max < 2.8 a small suppression can be 
seen. Final-state yield distributes in the central rapidity re- 
gion dominantly, thus production in central rapidity is af- 
fected by nuclear matter effects most intensively. 

Fixing the transverse energy of a jet in final state Eji = 
100 GeV, we survey the nuclear modifications for transverse 
momentum spectra of the other jet as displayed in Figure 
10. The yields in p+Pb and Pb+Pb collisions with different 
nPDFs are both enhanced, except that at very large Ej re- 
sults with EPS give a small suppression in Pb+Pb at ^/s = 
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Fig. 10 Nuclear modification factors for transverse energy spectra of 
dijet involving a fixed energy £"7-1 = 100 GeV , in the set of different 
nPDFs. 



Fig. 11 Nuclear modification factors of dijet asymmetry distribution 
for p+Pb collisions (top) and for Pb+Pb collisions (bottom) with dif- 
ferent nPDFs. 



2.76 TeV, which are in accordance with the modifications for 
dijet invariant mass spectra. The CNM effects increase the 
dijet yield in p+Pb collisions about 10% at high Ej and the 
yield in Pb+Pb collisions about20% at low Ej- Furthermore, 
we can see the different modifications brought by different 
sets of nPDFs obviously. The difference becomes slightly 
larger in Pb+Pb collision. 

From transverse momentum spectra of dijet , we can de- 
rive the nuclear modifications of dijet asymmetry A j due to 
CNM effects with the cuts Eji > 100 GeV for leading jet 
and Et2 > 25 GeV for subleading jet shown in Figure 11. 



We find that the nuclear modification on dijet asymmetry is 
similar with that on dijet angular distribution. The modifi- 
cation factors are close to unity and nearly independent of 
Aj over the whole range. The fluctuation at large Aj can be 
ignored since the yield at large Ay is very small as revealed 
in Figure 3. The dependence of modification effects on Ay is 
canceled by the fraction of transverse energy as in Eq. (7). 
The CNM effects have little influence on dijet asymmetry. 
It shows that the rather flat Ay distribution of dijet asym- 
metry measured by ATLAS and CMS in Pb+Pb collisions 
should be attributed to the final-state QGP effect. Thus the 
dijet asymmetry is an excellent signal to inspect the final- 
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State QGP effects as well as dijet angular distribution due to 
their insensitivity to the CNM effects. 



4 Conclusion 

In high-energy nucleus-nucleus collisions particles produc- 
tions at large transverse momentum have been extended from 
leading hadrons or vector bosons to reconstructed jets. In- 
teresting observations of dijet asymmetry in Pb+Pb by AT- 
LAS and CMS give a first successful example of jet mea- 
surements in heavy-ion collisions and may demonstrate the 
fingerprint of jet quenching in nuclear reactions at LHC for 
the first time. To put the study of dijet production in heavy- 
ion collisions on a more rehable base the investigation of 
CNM effects is indispensable. 

In the paper we focus on CNM effects on the dijet pro- 
ductions in relativistic heavy-ion collisions at LHC with the 
next-to-leading order pQCD by incorporating different nPDFs 
set by EPS, EKS, HKN and DS. The dijet angular distri- 
bution dapA(AA-)/apj^^AA)dX, dijet invariant mass spectrum 
d'^<^pA(AA) I dMjjd\y\,nax, final transverse energy spectrum of 
dijet d'^OpA(AA)/dETidET2 as well as dijet momentum im- 
balance dOpA^AA^ /<ypA(AA)dA J at LHC are calculated and the 
corresponding nuclear modification factors due to CNM ef- 
fects are given. We found that dijet angular distributions 
and asymmetry distributions are insensitive to the initial- 
state CNM effects and thus make themselves the excellent 
jet observables to investigate the final-state hot QGP effects 
such as jet quenching, while the dijet invariant mass spec- 
tra and transverse energy spectra show a conspicuous de- 
pendence on initial-state CNM effects. The overall enhance- 
ment of d^apA(AA)ldMjjd\y\max over a wide region of Mjj 
and d^ cr^^ ^^A) /dEj i dETi over a large region of Ej for p+Pb 
and Pb+Pb at LHC due to CNM effects is demonstrated, 
which is opposite to the suppressions of jet spectra at high 
Er due to initial-state CNM effects at RHIC and the final- 
state jet quenching effects both at RHIC and LHC. Visible 
difference between utihzing of nPDFs of EPS, EKS, HKN 
and DS is observed in the nuclear modification factors for 
dijet invariant mass and transverse momentum spectra in 
heavy-ion collisions, especially for large colliding system 
such as Pb+Pb, which makes d'^apA[AA)/dMjjd\y\max and 
d'^<^pA(AA) /dET\dET2 potential observables to distinguish dif- 
ferent parametrizations of nPDFs. The sensitivity of the nu- 
clear modification of invariant mass spectrum to jet radius 
R and rapidity y is also explored. It is illustrated that though 
dependence of 7?*^" on the jet size R is negligible, the varia- 
tion of R'^JJ with different rapidity regimes is relatively large 
in Pb+Pb colUsions. 
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